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SUMMARY 


We  obtain  for  several  geometric  situations  the  expression  for 
the  probability  of  insufficiently  accurate  data  on  the  velocity  of  a 
surface  target  as  a  function  of  the  parameters  scan  rate,  nearest 
approach  distance  of  an  observatory  satellite,  the  detection  radii  of 
the  radar  and  the  probability  of  detection  per  scan.  The  actual 
number  of  observations  taken  per  pass  of  the  satellite  is  random  and 
in  one  instance  we  give  explicitly  the  expectation  of  the  probability 
of  accurate  data.  This  expectation  can  be  calculated  in  other  cases. 

This  paper  allows  the  determination  of  the  capability  of  the 
radar  necessary  to  produce  an  acceptable  probability  of  obtaining 
accurate  velocity  information. 


1.  INTPODL'CTION 


Consider  an  observatory  satellite  moving  in  a  straight  line  across 
the  plane  which  at  any  given  time  observes  via  radar  a  sector  of  an 
annulus  on  the  surface.  This  annulus,  with  satellite  position  at  the 
center,  we  call  the  surveillance  ring.  This  ring,  we  consider,  encounters 
a  target  which  is  also  moving  linearly  in  the  plane  at  a  velocity  much 
less  than  that  of  the  satellite. 

It  is  easily  seen  that  if  we  regard  the  satellite  as  being  fixed,  the 
target  still  travels  linearly  across  the  surveillance  ring.  We  assume 
for  this  argument  that  the  satellite  is  fixed.  Note  that  from  this 
assumption  it  would  follow  that  the  surface  was  moving  which  would  have  to 
be  taken  into  account,  if  we  were  interested  in  estimating  the  satellite’s 
position  on  the  surface.  This  particular  problem  already  has  been  the 
subject  of  two  previous  documents,  [2]  and  (3),  and  we  are  now  interested  in 
a  different  one  namely,  the  probability  of  obtaining  good  data  from  the  radar. 

Before  we  proceed  further,  let  us  introduce  the  following  notation: 

Rj  is  the  minimum  detection  radius  of  the  radar, 
is  the  maximum  detection  radius, 

D  is  the  nearest  approach  distance  of  the  target  to  the  »'adar, 

p  is  the  probability  of  detection  per  scan  of  the  target  in 
the  surveillance  ring, 

h  is  the  distance  the  target  moves  relative  to  the  satellite 
between  scans. 
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'1'  is  tlie  distance  between  observations  resulting  in 

detection  that  must  be  obtained  to  procure  the  velocity 
data  needed, 

U  is  tile  distance  the  target  moves  into  the  annulus  before 
being  scanned. 

At  each  sweep  of  the  radar  across  the  target,  when  the  target  is  at 
least  but  less  than  a  distance  R2  away,  the  event's  detection  or 

non-detection  can  occur  with  probabilities  p  and  q  »  1  -  p, 
respectively.  Moreover,  successive  sweeps  of  the  radar  are  stochastically 
Independent  and  so  we  have  a  sequence  of  Bernoulli  trials  as  the  target 
moves  across  the  surveillance  ring.  We  assume  the  probability  of 
detection  is  zero  when  the  target  is  outside  the  surveillance  ring. 

2.  ONE-STAGE  ACCEPTANCE  WHEN  THE  DISTANCE  OF  NEAREST  APPROACH  EXCEEDS  THE 
MINIMUM  DETECTION  RADIUS 

We  now  assume  that  unsatisfactory  data  will  be  obtained  if  less  than 
two  detections  are  obtained  or  two  or  more  detections  are  obtained  but 
the  first  and  last  are  less  than  a  distance  T  apart. 

Whenever  Rj^  <  D  <  R^, 


we  have  the  following  situation: 
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The  length  of  the  ;ia:h  of  the  target  when  it  is  within  the 
surveillance  annulus  is  2d  and 


d  = 


The  number  of  scans  of  the  radar  is  N  where 


(2.0.8) 


(2.0.9) 
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■Aiiii  lx|  denoting  the  largest  Integer  less  than  x  for  x  any 
real  numher.  Because  we  have  a  sequence  of  N  Bernoulli  trials  each 
with  probability  of  success  p,  we  introduce  the  notation  for  the 
cumulative  binomial  probability 


B(k:n,p) 


0) 


pj(l  -  p)"-^ 


this  being  standard  since  the  classic  text  of  Feller  [1].  The 
probability  of  obtaining  less  than  two  detections,  that  is  either  zero 
or  one  detection,  is 


B(l:N,p)  =  +  Npq^ 


(2.1) 


Now  the  probability  that  we  have  at  least  two  detections  when  the  first 
detection  occurs  at  the  r^^  scan  and  the  last  at  the  r  +  s*”''  scan 
is  seen  to  be 

^  2  N-s-1 

TT‘’i"rPr+s  rT‘>i=P‘> 

1=1  i=r+s+l 

where  is  the  probability  of  detection  at  the  i^^  scan.  Hence, 

the  probability  of  at  least  two  detections  with  the  first  and  the  last 
less  than  or  equal  k  units  apart,  where 


k  =  [T/h], 


(2.1.1) 
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IS 


•^1  = 


k 

2  N-s-l 


.2  N-s-l 
,(N  -  s)p  q 


Letting  the  derivative  with  respect  to  the  variable  q  be 
denoted  by  D  , 

q 


by  using  the  fact  that  for  anv  positive  integers  N,k 


Now 


(2.1.5) 


q  -»•  (N-k)pq _ -  (q  -f-  Nq  ) 


Hence 


B(l:N-k.p)  -  B(l:N.p) 

2 

P 


(2.1.6) 


Qj  =  B(l:N-k,p)  -  B(l:N,p). 


(2.2) 


Thus  the  probability  ot  bad  or  insufficiently  accurate  data  for  the 
first  case  is  obtained  bv  adding  Equations  (2.1)  and  (2.2)  and  is 


=  B(l:N-k,p)  =  (N-k)pq^  ^  +  q^  ^ 


(2.3) 


where  N  is  given  in  (2.0.9),  k  is  given  in  (2.1.1). 
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W'v  llu’  distance  the  target  penetrates  into  the  surveillance  ring 
before  it  is  scanned  by  the  radar  is  a  random  variable  which  we  called 
U.  Thus,  any  function  of  U  is  also  random  e.g.  N  and  expressions 
and  Q2  above. 

If  we  assume  U  is  uniform  on  (0,h)  and  we  let 
[d/h]  =  n  (d/h)-n  =  f 

be  the  integral  and  fractional  part  of  d/h,  one  can  see 


N 


n  +  1 

n 


with  probability  f 

"  "  1  -  f. 


(2.3.5) 


Thus,  we  have  by  taking  the  expectation,  the  expected  probability 
of  obtaining  unsatisfactory  data 


EQ2  =  fB(l:n-k+l,p)  +  (1-f )B(l:n-k,p) .  (2.3.6) 

Now  this  form  lends  itself  to  easy  use  with  tables  of  the  cumulative 
binomial  distribution  e.g.  the  Harvard  tables  [4]. 

One  can  obtain  alternate  forms  by  manipulation  e.g. 

EQ2  =  fp(n-k)q"~*^  +  (1-f)  (n-k)pq"  ^  +  q^‘ 


and  again 

EQ2  •  ((n-k)p(fq  +  1  -  f)  +  q)q"  ^  ^  (2.4) 

as  a  final  form  for  the  probability  sought  for  Case  I. 
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3. 


acceptance  utien  the  nearest  approach  distance  is 

THE  MINIMl'M  DETECTION  RADII’S 


LESS  THAN 


We  now  consider  the  case  when  0  <  D  ^  R^.  At  this  distance,  the 
relative  motion  of  the  target  to  the  satellite  carries  it  inside  the  blind 
range  of  tlie  radar.  The  situation  is  as  follows: 


MOTION  OF  TARGET  RELATIVE  TO  SATELLITE 


We  now  number  the  scans  in  the  order  of  their  occurrence  inside  the 
surveillance  ring  from  the  top  down.  And  we  write 

Pj  =  probability  of  detection  at  the  scan  j  =  l,...,N 


where 
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^2  =  ‘"^2  ■ 


-M*' 

J,  =  /r?  - 


(3.1) 


Also  let  N’  be  the  number  of  observations  obtained  in  passing  through 
the  first  segment  of  the  ring. 

=  [  (d2  ~  -  U)/h]  +  1  (3.2) 

and  N  +  1  be  the  number  of  the  first  observation  obtained  in  the 
passage  through  the  second  segment  of  the  ring. 


N"  =  [(dj  +  d2  -  U)/h]  +  1. 
Since  in  this  case  we  have 


-i: 


j=N'+l,...,N" 
P  otherwise 

The  probability  of  zero  detections  is 


N 


IT 

i=l 


N'+N-N" 


The  probability  of  one  detection  is 


i  Pj  TT  ‘’i  =  ^  ^  pqN'+N-N"-l 

“  jfr+1 


(3.3) 


(N*  +  N  -  N")pqf'  +N-N--1  _ 


(3.6) 
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Hence,  as  before  we  see  the  probability  of  less  than  two  detections  is 


B(1:N'  +  N  -  N",p). 


(3.6.5) 


The  probability  of  two  or  more  detections  with  the  first  detection 
at  tiie  r*^^  scan  and  the  last  detection  at  the  r  +  s*^^  scan  is 


r-1 


9  •  P  P  II  9 .  = 

.  ,  1  r'^r+s  .  '  '  ^i 

1=1  i=r+s+l 


2  N-s-N"+N’-l 


P  9 


0 


for  r+s  ^  N' 
or  r  ^  N"+l 
otherwise . 


Hence,  the  probability  of  at  least  two  detections  with  the  first  and 
last  detections  less  than  or  equal  k  units  apart,  with  k  =  [T/h],  is 


‘’2'^  TT  ’i)  • 

s=i  /  yi=r+s+l  / 


-ss 


^2  = 


2  N-s-N"+N-1 
P  9 


(|r  ^  N"+lJ  +  |r  +  s  N'l) 


where  fx  ^  t|  is  the  indicator  function  of  the  relation  taking  the 
value  one  for  all  values  of  the  variable  for  which  it  is  true  and  the 
value  zero  otherwise.  Now  we  have 


„i  X  ■*  2  N-s-N"+N'-1 
^2  =  >  P  9 


((N-s-N")’’  +  (N'-s)'^ 


where  again  we  have  introduced  the  notation 


!x  if  X  >  0 

0  if  X  <  0 


(x) 


U) 


2  N 

=  p  f| 


g 


(K-s-N")  q 


N-s-X"-l 


2  N-N"' 
P  f] 


s=i 


(N'-s)  q 


N'-s-l 


where  in  this  paper  we  make  the  convention  that 

{a,b}  =  r.ln(a,b) . 

Now  using  the  result  (2.1.6)  we  have 


(3.7) 


0^  =  q^' (B(l:N-N’’-{k,N-N"}  ,p)  -  B(1:N-N",p)) 

+  q^"^"(B(l:N-{k,N’},p)  -  B(l:N'.p)) 

and  by  adding  (3.6.5)  and  G2  above  we  have  the  probability  of  bad 
data  for  this  case,  namely, 

=  B(1:N+N'-N",p)  +  q^'(B(l:N-N"-{k,N-N"},p)  -  B(1:N-N",p)) 

+  q^"^"(B(l:N'-{k,N'},p)  -  B(l:N’-p)). 


But  we  notice  that 

B(1:N+N'-N",p)  -  q^'B(l:N-N’’,p)  -  q^"^"B(l :N ’ , p) 
N-N"+N ' 

=  -q 


(3.7.1) 


hence,  the  answer  is 
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•^3  = 


N' 

q  B(l:N-N"-{k,N-N"},p)  + 


N-\'"  V' 

q'"  B(l:N'-{k.,N'},p)-q'^  ' 


N"+N 


(3.8) 


To  find  the  expected  value  of  this  expression,  we  are  reduced 
to  a  tedious  evaluation  of  special  cases.  Take  as  the  integral  and 
fractional  parts  of  the  following: 


dj/h  =  n.  +  f^  i=l,2, 

with  V  =  U/h  being  a  random  variable  uniform  on  (0,1), 


1  2n2 

if 

0 

>  2f2  - 

V  >  -1 

N  =  [2n2+2f2-V] 

+  1  = 

1  ^"2 

+  1  if 

1 

>  2f2  -  V 

>  0 

1  2n2 

+  2  if 

2 

>  2f2  -  V 

>  1 

"2  ■ 

+  1 

if 

l>f2  - 

fl  -  V 

>  0 

[n2-n^+f2-f  j^-V] 

+  1  = 

"2  - 

"1 

if 

0>f2  - 

^  -  V 

>  -1 

1  ‘^2 

"1  -  ^ 

if 

-1  >  f2 

-  ^1  - 

V  > 

n  +  n  +  2  if  2  >  f,  +  f^  -  V  >  1 


N"  =  [n^+n2+fj+f2-V]  +1=  |n^+n2+l 


"l  "2 


if 


1  >  f^  +  f2  -  V  >  0 


if  0  >  f^  +  f^  -  V  >  -1 


Now  the  evaluation  of  the  expectation  of  (3.8)  with  respect  to  the 


random  variables  N,N',N" 


would  be  so  tedious  to  write  out  and  would 
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coiivov  such  little  information  additional  to  (3.8) 
from  doinj;  so. 


that  we  refrain 


4.  THE  TWO-STAC.E  ACCEPTANCE  WITH  NEAREST  APPROACH  EXCEEDING  THE  MINIHUM 
DETECTION  RADIUS 

We  adopt  the  same  probabilistic  model  as  before  only  now  the  criterion 
for  insufficient  data  we  take  as  being  less  than  two  detections  or  two 
detections  separated  by  at  most  a  distance  T^  or  three  or  more  detections 

separated  by  at  most  T^. 

As  before,  we  have  in  the  first  case  the  probability  of  less  than 
two  detections  being 


B(l:N,p) 


+  Npq 


N-1 


(4.1) 


The  probability  of  exactly  two  detections  with  the  first  and  last 
less  than  or  equal  =  [Tj^/h]  units  apart  is 


k,  N-s 


2  N-2 
p  q 


(N-s)p  q 


2  N-2  k^(2N-k^  1)  2  jj-2 


p  q 


(4.2) 


s=l  r= 


The  probability  of  three  or  more  detections  with  the  first  and 
last  detections  less  than  or  equal  k2  =  (T2/h]  units  apart  is  easily 
seen  to  be  by  comparison  with  the  argument  of  Section  2 


2 

P 


(1-q 


1.  N-s-1 

)q 


k 

2xr^  .  .  N-s-l  N-2 
p  y  (N-s)(q  -q 

s= 
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N-s\  ^^2  2  N-2 

I  -  i  P  ^ 

(k  -l)(2N-k  -2) 

=  B(l:N-k2,p)  -  B(l:N-l,p)  -  p  q  — - - - ^ -  .  (4.3) 

Thus,  the  probability  of  insufficient  data  is  the  sum  of  (4.1),  (4.2) 
and  (4.3)  which  is 


k  (2N-k  -1)  2 

=  B(l:N,p)  +— - -  p^q‘  +  B(l:N-k2,p)  -  B(l:N-l,p) 

-  P  q  - ^ - 

since  B(l:N,p)  -  B(l:N-l,p)  =  -(N-l)p^q^  ^ 

(4.4) 

By  letting  k^  =  k.^,  we  see  this  reduces  to  the  preceding  result 
of  Equation  (2.3).  We  do  not  bother  with  the  expert  it  ion  of  f',  . 
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C.  = 

1 


{k^,N' }(2N'-{k^,N' }-l)  +  {k^,N-N"}(2N-2N"-{k^,N-N"}-l) 


,  _  2  N'+N-N”-2 

Q4  -  q 


The  probability  of  three  or  more  detections  with  the  first  at  the 
r^^  scan  and  the  last  at  r  +  s^^  scan  by  analogy  with  previous  work  is 


2,,  s-1.  N'+N-N"-s-l 

P  (1-q  )q 


I  0 


if  r  +  s  N ' 
or 

r  2  N"  +  1 
otherwise. 


Therefore,  the  probability  of  three  or  more  detections  with  the  first 
and  last  less  than  or  equal  k^  units  apart  is 


k^  N-s 


2,,  s-1,  N'+N-N”-s-l 

p  (1-q  )q 


(fr  +  s  1  N'l  +  tr  2  N”  +  1^ 


s»l  r»l 


^2 

S=1 


((N'-s)^  +  (N-s-N")  ) 


'+N-N"-s-l  /  N'+N-N"-2  N'+N-N"-s-l 


.  "  ‘-(N'-s)  q .  V(N-N"-s)  q‘ 


s=  1 


VH  x'  N'+N-N"-2\ 
-(N-N  -s)  q  } 
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=  P 


{k^.N-N"} 

.  pV'd  (  2  ^'‘-'‘''-=)- 

S=1  '  '  S*1  ' 


„  2  N'+N-N"-2 

^>2?  q 


q"'"  (B(l:NMk2.K'),p)-B(l:N',p))-C2P^'‘'+''-'‘''-2 

+  q"  (B(l:N-K''-(k2,»-ll'')  ,p)-B(liN-N",p)). 


Hence,  the  probability  that  we  have  insufficiently  accurate  data  for 
this  case  is 


Q5  =  B(l:N'+N-N".p)  +  p2qN*+N-N"-2 

+  q”"”  (B(l:Ni{k2.N'} ,p)-B(l:N',p)) 
N'  / 

+  q  (B(l;N-N"-{k  ,N-N"},p)-B(1:N-N" 


(C1-C2) 


which  by  the  simplification  of  (3.7.1)  is 


writing  out  the  expectation  of  Q^. 
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